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Abstract 

We obtain the generating function for shot noise through a disordered multi- 
channel conductor in the zero temperature quantum regime. The distribution of 
charge transmitted over a fixed time interval is found to be approximately Gaussian. 
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The nature of shot noise in mesoscopic systems at low temperatures has attracted 
considerable interest in recent years. Lesovik [Q] showed that the intensity of the shot 
noise in a single channel two terminal conductor is supressed in the quantum regime 
compared to its classical value. This result has been generalized and applied to various 
other cases || [3|, §]]. Levitov and Lesovik J5J later was able to obtain the entire distribution 
function for the charge transport over a fixed time interval for a single channel conductor 
and showed that the distribution function becomes binomial in the quantum regime as 
opposed to Poissonian in the classical regime. This result has been generalized to include 
shot noise in a superconductor-normal metal point contact as well. 

In the present work we generalize the work of ref. || and evaluate the shot noise 
distribution function for a multichannel conductor. We obtain the generating function for 
the distribution of charge transmitted over a time interval t and show that the charge 
distribution is approximately Gaussian. We find that in some limit the shot noise is 
reduced by a factor close to 1/3 compared to the Poisson value as obtained in ref. ||, 
but we show that there is a correction term independent of the value of conductance that 
can be important in the large t limit. Our results disagree with a recent calculation by 
Lee et al |7]], and we comment on the differences. 

When a dc voltage V is applied across the system in a two terminal geometry, the shot 
noise generating function for a multichannel conductor may be written as a product over 
generating function for independent channels, each of which has a binomial distribution 

& 

X (A)=n(^ + l-T a ) M ; z = e l \ (1) 

a 

where M = (e/h)Vt can be interpreted as the number of attempts in time to, and T a is 
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the transmission probability corresponding to channel a. The probability distribution for 
Q charges [each of charge e] transferred in a time to is 

P(Q,t Q )= d\e* QX < X (\,t )> (2) 

J —IT 

where the angular bracket represents an average over the distribution of the T a 's in disor- 
dered mesoscopic conductors. One can also read off the moments or the cumulants of the 
charge distribution directly by expanding < x(A,t ) > or In < ^(A, t ) > in power series 
of A. 

It has been shown || that the distribution of the transmission probability in the diffu- 
sive regime for an iV-channel disordered conductor can be obtained from an appropriate 
random matrix theory constructed for the matrix X = [TT^ + (TT^)^ 1 — 21] /A where T is 
the 2N x 2N transfer matrix describing the conductor and / is the identity matrix. The 
N doubly degenerate real non-negative eigenvalues x a of the matrix X are related to the 
transmission probability T a by the relation 

T a = — - — ; < x a < oo. (3) 

1 + Xa 

According to the random matrix theory 0, the joint probability distribution for the N 
eigenvalues then be written in the form 

p(xi---x N ) = \{\x a -x b \ a \{e- v ^\ (4) 

a<b c 



where V(x) is a Lagrange multiplier function |T(J that incorporates physical constraints 



like a given value of the first moment, and a is a symmetry parameter. The value of a 
is 1,2 or 4 depending on whether the symmetry of the matrix ensemble is orthogonal, 
unitary or symplectic. A reasonably good description of a disordered conductor in the 



diffusive regime M, [11], is obtained from the choice V{x) = tx, where the parameter 
t is related to the dimensionless ohmic conductance go of the system by t = ag\j1N ' . 
The conductance g is given by the relation g = J2 a T+^~; an d is therefore a linear statistic 
of the eigenvalues x a ; the variance of such linear statistics have well known properties 
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The generating function for shot noise is, however, not a linear statistic; it is a product 
statistic as given by eq. (1). To our knowledge, fluctuation properties of such product 
statistics within a random matrix framework have not been calculated before. [One might 
attempt to avoid the problem by evaluating the average < ln% >, which is a linear statis- 
tic J7|, instead of the appropriate In < \ >; this is a reasonable approximation only if 
the distribution is very sharply peaked. Our results differ from those of M]. We will 
evaluate the average of the product statistic within the random matrix framework using 
a continuum approach [ |T5[] , which is known to work very well in the diffusive regime. In 
this approach the probability distribution (4) is written as the exponential of a fictitious 
Hamiltonian where the eigenvalues are considered as particles repelling each other loga- 
rithmically as in a Coulomb fluid, within a confining potential V(x), at a temperature 
1/a. 

In terms of the variable x a , the generating function can be rewritten as 

x (X,t ) = e w E.M(*.+*)/(*.+i)]. ( 5 ) 

The average of x over the distribution of then be written as 

(X(A)> = Zm/Zo = e-( F -- F °) (6) 
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where 



Jl dx a exp 



a^ln|x a - x b \ - tJ2 x a + M^ln 

a<b 



X a -\- Z 



a 



(7) 



X a + 1 

Here Z , defined as Zm=o is the partition function and F Q , defined as Fm=o is the free 
energy for the coulomb fluid characterized by the logarithmic repulsion and a linear con- 
fining potential as shown in eq. (7) for M = 0. It is well known that for a linear confining 
potential the density of eigenvalues at the origin scales with the number of eigenvalues N. 



Therefore the large N continuum approximation is valid [fLql , and we can use the coulomb 
fluid results 



U' M (x) - a 



dy 



'i x — y 

where &m{x) is the 'effective' density of eigenvalues contained in an interval J, the prime 
denotes a derivative with respect to x, and 



U M (x) = U {x) + U c {x) =V(x)-Mhi 



X + z 



x + 1 



(9) 



is the 'effective' confining potential, which includes the confining potential Uq{x) = 
Um=o(x) = V(x), as well as the 'correction' U c (x) = Mln 2±| . [We use the word 



'effective' within a quotation mark because the parameter z is complex, and the density 
or the confining potential has only a mathematical meaning.] The range of integration I is 
determined by the allowed range of eigenvalues as well as the normalization requirement. 



In terms of the density, the free energy Fm can be obtained [16| from the relation 



F M = (1/2) A M N + J dxa M (x)U M (x), 



where Am is the chemical potential determined from the integration constant, 



(10) 



U M (x) - a J dya M (y) In \x - y\ = A M . 



11) 



The general solution to equation (8) for the density is given by [I = (0, b)} 



*-<*> = ii/^ /'A&w, (12) 



7r 2 a V a; Jo y — xVb — y 

where b is the upper limit for the density set by the symmetry and normalization require- 
ments. 

To obtain explicit expressions, we write 

<?m{x) = vo(x) + <?c{x), (13) 

where cxo(x) = <Jm=o(x) is the density for the potential Uq(x) alone, and a c (x) is the cor- 
rection to the density due to U c (x). Then for Uo(x) = tx, and for non-negative eigenvalues, 
we obtain 



*o(*) = — J— > ( 14 ) 
na V x 



and 



, M(l-v^) x-Vi . . , 1E . 

^ x = ^—7 — nr? — ar9Z < n - 15 

The normalization condition J 6 cIx<jm(x) = N gives us the value for b, iV = tb/2a. Note 



that the contribution from a c to the normalization integral is zero [17], 



oo 







dxa c (x) = 0. (16) 



Therefore we may conclude that the chemical potential Am differs from Aq by a negligible 
amount as iV — > oo. This provides a major simplification, and the change in free energy 
is then given by 

F M -F = (1/2) / dxa U c + (1/2) / dxa c U + (1/2) / dxa c U c , (17) 

JO JO JO 
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where oo, Uo, cr c and U c are as denned before. In the limit of large system size, b = 
4N 2 /g 2 — > oo (go = Nl/L, where / is the mean free path and L is the system size), all 
integrals can be explicitly evaluated, and the generating function is given by 



(xO)) = exp 



-Mg (l -^) + ^(l- z ) + —\n 



, z = e iX , |A| < 7r 
(18) 



where we have used the conductance g instead of the parameter t. This is our principal 
result. Again for large system size, << Mg , and the term proportional to (1 — z) can 
be neglected compared to the term proportional to (1 — y/z). However, for long enough 
measurement time to (a single energy generating function can be considered only in this 
limit), the term proportional to M 2 may not be negligible (M = eVt /h). Note that this 
term is independent of the conductance go, but depends on the symmetry parameter a. 

We can use eq. (2) to evaluate an approximate probability distribution P(Q,t ). We 
use the fact that major contributions to the integral in (2) come from regions around 
A = where expansions in powers of A can be used. If we keep terms upto order A 2 , the 
limits of the integral can then be extended to ±oo. The result, valid around the mean 
Qo, is a Gaussian distribution: 

(Q ~ Qo) 2 



P(Q,t ) ~ exp 
where the mean value 



2(* + £) 



(19) 



< Q >= Qo = ^^ M 9o- (20) 
We can also expand In (x) in power series of A to obtain directly the cumulants of the 
charge distribution. For — > 0, we get 

i/\ n M (Qo M M 2 \x 2 (Q \x* (Qo M 2 U 4 
fo{x)=Qox+\— + — )^+[ — )^+\ — - — \-7 + ... (21) 



2 8a J 2! V 4 / 3! \ 8 64a/ 4! 



where the cumulant << Q k » can simply be read off as the coefficient of x k /k\ in the 
expansion (21). The first two cumulants agree with the approximate Gaussian distribution 
obtained above. Note that Qo defines the average charge transferred in time to, while the 
intensity of the shot noise is given by << Q 2 ». A completely uncorrelated classical 
system would have a Poisson distribution |§ with < Q > as well as << Q 2 >> given 
by Qp = Mg . If the term proportional to M 2 is neglected, the shot noise is reduced 
compared to Qp by a factor which is close to 1/3 obtained in ref. ||. However 
the correction term, which is independent of go, may become important for long enough 
observation time and/or stronger disorder. On the other hand, for stronger disorder, the 
appropriate random matrix model [II] has a confining potential V(x) ~ [lnx] 2 , for very 



large x, and the density near the origin no longer scales with N. In this regime the usual 
large N continuum approximation breaks down. The shot noise distribution in this regime 
is under investigation. 
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